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Abstract  We  consider  the  broadcasting  problem  in  multi -radio  multi-channel  ad  hoc 
networks.  The  objective  is  to  minimize  the  total  cost  of  the  network-wide  broadcast, 
where  the  cost  can  be  of  any  form  that  is  summable  over  all  the  transmissions  (e.g., 
the  transmission  and  reception  energy,  the  price  for  accessing  a  specific  channel). 
Our  technical  approach  is  based  on  a  simplicial  complex  model  that  allows  us  to  cap¬ 
ture  the  broadcast  nature  of  the  wireless  medium  and  the  heterogeneity  across  radios 
and  channels.  Specifically,  we  show  that  broadcasting  in  multi-radio  multi-channel  ad 
hoc  networks  can  be  formulated  as  a  minimum  spanning  problem  in  simplicial  com¬ 
plexes.  We  establish  the  NP-completeness  of  the  minimum  spanning  problem  and 
propose  two  approximation  algorithms  with  order-optimal  performance  guarantee. 
The  first  approximation  algorithm  converts  the  minimum  spanning  problem  in  sim- 
plical  complexes  to  a  minimum  connected  set  cover  problem.  The  second  algorithm 
converts  it  to  a  node-weighted  Steiner  tree  problem  under  the  classic  graph  model. 
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These  two  algorithms  offer  tradeoffs  between  performance  and  time-complexity.  In 
a  broader  context,  this  work  appears  to  be  the  first  that  studies  the  minimum  span¬ 
ning  problem  in  simplicial  complexes  and  weighted  minimum  connected  set  cover 
problem. 

Keywords  Broadcast  •  multi-radio  multi-channel  •  ad  hoc  network  •  simplicial 
complex  •  minimum  spanning  •  minimum  connected  set  cover. 


1  Introduction 

Multi-Radio  Multi-Channel  (MR-MC)  wireless  networking  arises  in  the  context  of 
wireless  mesh  networks,  dynamic  spectrum  access  via  cognitive  radio,  and  next- 
generation  cellular  networks  [1 1].  By  the  use  of  multiple  channels,  spatially  adjacent 
transmissions  can  be  carried  over  non-overlapping  channels  to  avoid  mutual  interfer¬ 
ence.  Furthermore,  each  node,  equipped  with  multiple  radios,  is  capable  of  working 
in  a  full-duplex  mode  by  tuning  the  transmitting  and  receiving  radios  to  two  non¬ 
overlapping  channels. 

The  increasing  demand  for  high  data  rate  and  the  persistent  reduction  in  radio 
costs  have  greatly  stimulated  research  on  MR-MC  networks.  Considerable  work  has 
been  done  on  capacity  analysis,  channel  and  radio  assignment  [10,  16,  18,  21],  and 
routing  protocols  [10,  18].  In  this  paper,  we  consider  the  broadcasting  problem  in 
MR-MC  ad  hoc  networks. 


1.1  Broadcasting  in  Single-Radio  Single-Channel  Networks 

Broadcasting  is  a  basic  operation  in  wireless  networks  for  disseminating  a  message 
containing,  for  example,  situation  awareness  data  and  routing  control  information,  to 
all  nodes.  For  a  Single-Radio  Single-Channel  (SR-SC)  network,  a  key  question  for 
the  network-wide  broadcast  is:  which  set  of  nodes  should  be  selected  to  transmit  such 
that  the  total  cost  (such  as  energy  consumption  or  the  number  of  transmissions)  is 
minimized.  In  contrast  to  the  wireless  broadcast  problems  for  minimizing  the  energy 
consumption  and  the  number  of  transmissions  which  are  shown  to  be  NP-complete 
in  [22],  their  counterparts  in  wired  networks  have  polynomial  solutions. 

The  complexity  of  the  problem  arises  from  the  broadcast  nature  of  the  wireless 
medium:  a  single  transmission  from  one  node  can  reach  all  the  other  nodes  within  the 
transmission  range  of  this  node,  but  it  may  cause  interference  to  other  nearby  trans¬ 
missions.  This  ''node-centric"  nature  of  the  wireless  broadcasting  problem  along  with 
the  mutual  interference  between  concurrent  transmissions  complicates  the  design  of 
efficient  broadcasting  algorithms. 


1.2  Broadcasting  in  MR-MC  Networks 

In  an  MR-MC  ad  hoc  network,  such  as  the  DARPA  Wireless  Network  after  Next  (W- 
NaN)  [19],  each  node  is  equipped  with  multiple  radios,  each  operating  on  a  different 
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channel.  The  introduction  of  multiple  channels  and  multiple  radios  further  compli¬ 
cates  the  design  of  an  efficient  broadcasting  scheme.  Since  the  number  of  radios  at 
a  node  is  usually  smaller  than  the  number  of  channels,  the  broadcast  scheme  should 
decide  not  only  which  nodes  act  as  relays  but  also  for  those  relay  nodes,  which  chan¬ 
nels)  should  be  assigned  to  the  transmitting  radio(s).  Given  the  selection  of  the  relay 
nodes,  two  simple  broadcast  schemes  are:  (i)  transmitting  multiple  copies  of  the  mes¬ 
sage  on  all  channels;  (ii)  transmitting  a  single  copy  of  the  message  on  a  common 
channel  dedicated  to  broadcasting.  Both  schemes  are  inefficient.  For  the  latter  one, 
if  the  broadcast  load  is  high,  the  common  channel  will  be  overwhelmed,  even  while 
plenty  of  other  channels  are  free. 

One  subtle  issue  is  the  complication  of  the  wireless  broadcast  advantage.  In  an 
MR-MC  network,  if  the  radios  of  the  neighboring  nodes  are  tuned  to  different  chan¬ 
nels,  a  single  transmission  on  one  channel  cannot  reach  all  the  neighboring  nodes 
simultaneously.  In  other  words,  only  the  neighboring  nodes  on  the  same  channel  can 
share  the  wireless  broadcast  advantage.  More  precisely,  the  concept  of  neighborhood 
must  be  defined  both  by  radio  range  and  channel.  Another  subtle  issue  is  channel  het¬ 
erogeneity.  Channels  may  have  different  bandwidth,  fading  condition,  and  accessing 
cost,  leading  to  different  implications  for  the  total  broadcast  cost. 

Broadcasting  in  MR-MC  networks  is  thus  a  multi-faceted  problem,  involving 
channel  assignment,  relay  node  selection,  and  channel  selection  for  the  source  and 
relay  nodes.  In  this  paper,  we  focus  on  the  latter  two  issues  by  assuming  a  given 
channel-to-radio  assignment.  To  avoid  the  hidden  channel  problem  [16],  two  nodes 
that  are  two-hops  away  from  each  other  are  assigned  two  distinct  sets  of  channels. 
Our  design  objective  is  to  minimize  the  total  broadcast  cost,  where  the  cost  can  be 
of  any  form  that  is  summable  over  all  the  transmissions,  including,  for  example,  the 
transmission  and  reception  energy1,  the  price  for  accessing  each  channel. 


1.3  A  Simplicial  Complex  Model  for  Broadcasting  in  MR-MC  Networks 

Our  technical  approach  is  based  on  a  simplicial  complex  model  of  the  broadcasting 
problem  in  MR-MC  networks.  A  simplicial  complex  is  a  collection  of  nonempty 
sets  with  finite  size  that  is  closed  under  the  subset  operation.  In  other  words,  if  a 
set  s  belongs  to  the  collection,  all  subsets  of  s  also  belongs  to  the  collection.  An 
element  of  the  collection  is  called  a  simplex  or  face.  This  subset  constraint  is  often 
satisfied  in  the  network  context.  For  example,  subsets  of  a  broadcast/multicast  group 
are  broadcast/multicast  groups,  subsets  of  a  clique  are  cliques.  A  simple  example 
of  graph  and  simplicial  complex  is  given  in  Fig.  1.  While  the  concept  of  simplicial 
complex  has  been  around  since  the  1920’s,  many  well-solved  fundamental  problems 
in  graph  remain  largely  open  under  this  more  general  model. 

We  use  a  simplicial  complex  model  rather  than  a  graph  because  the  simplicial 
complex  more  naturally  captures  the  broadcast  channel,  and  the  distinction  and  dis¬ 
jointness  between  broadcasting  on  different  channels.  Further,  costs  can  be  attached 
to  faces  (simplices)  in  a  way  not  easily  possible  with  graphs. 


1  The  ‘reception  energy’  denotes  the  energy  consumed  by  the  radio  in  reception  mode. 
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Consider  an  example  MR-MC  network.  As  shown  in  Fig.  2,  after  the  channels 
are  assigned,  the  network  is  partitioned  into  cliques  of  nodes.  Such  clique  assignment 
scheme  is  implemented  in  the  DARPA  WNaN  system  and  was  used  in  an  exercise 
at  Ft.  Devens  [19].  A  clique  consists  of  the  nodes  which  are  within  each  other’s 
transmission  range  and  share  at  least  one  common  channel,  and  two  cliques  are  s- 
pliced  ( i.e connected)  via  nodes  operating  on  multiple  channels  shared  in  common 
by  the  two  cliques.  Within  each  clique,  depending  on  the  cost  function,  the  transmit¬ 
ter  decides  which  dimension  simplex  (i.e.,  a  subclique  or  the  clique  itself)  in  a  clique 
complex  to  activate.  The  message  for  the  network-wide  broadcast  is  thus  propagated 
through  a  sequence  of  cliques,  possibly  of  different  dimensions.  Note  that  the  unicast 
case  corresponds  to  a  clique  of  dimension  1  (an  edge).  This  example  could  also  ap¬ 
ply  to  the  case  where  nodes  may  have  multiple  radios,  perhaps  of  different  modality 
(e.g.,  RF  and  optical);  in  this  case,  there  may  also  be  a  cost  associated  with  switching 
modes. 

The  network-wide  broadcast  problem  can  be  formulated  as  the  minimum  span¬ 
ning  problem  in  simplicial  complexes.  A  clique  in  the  MR-MC  network  is  modeled 
as  a  simplex  in  the  simplicial  complex  (see  Fig.  2),  and  since  a  subset  of  a  clique 
is  still  a  clique,  the  constructed  simplicial  complex  meets  the  requirement  of  being 
closed  under  the  subset  operation.  The  minimum  spanning  problem  in  a  simplicial 
complex  is  to  find  a  connected  subset  of  simplices  that  covers  all  the  vertices  with 
the  minimum  total  weight,  i.e.,  the  Minimum  Connected  Spanning  Subcomplex  (M- 
CSSub)2.  Then  the  solution  to  the  network-wide  broadcast  problem  can  be  obtained 
by  solving  the  MCSSub  problem. 


1 .4  Minimum  Spanning  Problem  in  Simplicial  Complexes 

The  minimum  spanning  problem  in  a  graph  is  to  find  a  connected  subgraph  that 
covers  all  the  vertices  with  minimum  total  weight.  The  solution  must  be  a  tree  for 
graphs  with  nonnegative  weights  (hence  called  the  Minimum  Spanning  Tree  (MST)). 
There  are  several  polynomial-time  algorithms  for  MST,  e.g.,  Kruskal’s  Algorithm  and 
Prim’s  Algorithm  [9].  The  MST  problem  has  many  applications  in  network  planning, 
broadcasting  in  communication  networks,  touring  problems,  and  VLSI  design  [1]. 

With  the  addition  of  high  dimensional  simplices,  the  minimum  spanning  problem 
in  a  simplicial  complex  is  fundamentally  different  and  much  more  difficult  than  its 
counterpart  in  a  graph.  First,  unlike  the  case  in  a  graph,  the  MCSSub  of  a  simplicial 
complex  may  not  be  a  “tree”3.  As  illustrated  in  Fig.  3,  the  MCSSub  of  the  simplicial 
complex  is  the  three  filled  triangles  which  form  a  cycle.  Second,  while  simple  greedy- 
type  polynomial-time  algorithms  exist  for  finding  the  minimum  spanning  tree  in  a 


2  Strictly  speaking,  a  subcomplex  should  also  be  closed  under  the  subset  operation,  but  without  loss  of 
generality,  we  do  not  include  this  condition  in  the  definition  of  minimum  connected  spanning  subcomplex, 
which  is  also  more  relevant  to  the  broadcasting  problem  at  hand. 

3  Although  there  is  no  unified  definition  of  tree  in  simplicial  complexes,  a  couple  of  definitions  can  be 
obtained  by  generalizing  those  equivalent  definitions  of  tree  in  a  graph.  For  example,  simplicial  trees  can 
be  defined  based  on  the  universal  existence  of  leaves  in  any  subgraph,  or  the  uniqueness  of  simplicial  facet 
paths  (see  Sec.  2). 
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graph,  the  minimum  spanning  problem  in  a  simplicial  complex  is  NP-complete  as 
established  in  this  paper  (see  Sec.  3.1). 

We  develop  polynomial-time  approximation  algorithms  for  the  minimum  span¬ 
ning  problem  in  simplicial  complexes.  We  propose  two  algorithms:  one  reduces  this 
problem  to  a  minimum  connected  set  cover  problem,  and  the  other  reduces  the  prob¬ 
lem  to  a  node-weighted  Steiner  tree  problem  in  a  graph  derived  from  the  original 
simplicial  complex.  We  also  establish  the  approximation  ratios  of  the  two  algorithm- 
s.  Both  are  shown  to  be  order-optimal.  The  time-complexity  of  these  two  algorithms 
is  also  analyzed,  illustrating  the  tradeoff  between  performance  and  complexity  of¬ 
fered  by  these  two  algorithms.  In  a  broader  context,  this  work  appears  to  be  the  first 
that  studies  the  minimum  spanning  problem  in  simplicial  complexes  and  weighted 
minimum  connected  set  cover  problem. 


1 .5  Related  Work 

Broadcasting  in  MR-MC  networks,  mostly  in  the  context  of  wireless  mesh  networks, 
has  been  studied  for  different  optimization  objectives  (see  [2,  3,  12,  16]  and  ref¬ 
erences  therein).  Different  from  the  previous  work,  the  optimization  objective  in  our 
work  can  be  any  cost  function  which  is  summable  over  all  the  transmissions,  thus  tak¬ 
ing  into  account  channel  heterogeneity  (e.g.,  transmissions  on  different  channels  may 
consume  different  amounts  of  energy,  due  to  different  bandwidths  or  different  prop¬ 
agation  characteristics  or  some  other  factor).  We  point  out  that  neither  minimizing 
the  total  number  of  transmissions  nor  minimizing  the  total  number  of  radios  used  in 
the  broadcast  is,  in  general,  equivalent  to  minimizing  the  total  energy  consumption. 
The  reception  energy  is  ignored  if  the  total  number  of  transmissions  is  minimized, 
while  the  transmission  energy  and  the  reception  energy  are  equated  if  the  total  num¬ 
ber  of  radios  is  minimized.  More  importantly,  channel  heterogeneity  is  not  addressed 
if  these  two  objectives  are  optimized. 

Furthermore,  to  our  best  knowledge,  our  work  is  the  first  to  adopt  simplicial  com¬ 
plexes  to  model  and  solve  the  broadcast  problem  in  wireless  ad  hoc  networks.  For 
a  more  detailed  discussion  on  the  potential  applications  of  simplicial  complexes  in 
communication  and  social  networks,  readers  are  referred  to  [17]. 


2  Basic  Concepts  in  Simplicial  Complexes 

In  this  section,  we  introduce  several  basic  concepts  in  simplicial  complexes  [15]. 

An  (abstract)  simplicial  complex  is  a  collection  A  of  nonempty  sets  with  finite 
size  such  that  if  A  €  A,  then  V  B  C  A,  B  e  A,  i.e.,  A  is  closed  under  the  operation 
of  taking  subsets.  The  element  A  of  A  is  called  a  simplex  of  A;  its  dimension  (denoted 
by  dim  A)  is  one  less  than  the  number  of  its  elements.  Each  nonempty  subset  of  A 
is  called  a  face4  of  A.  The  dimension  of  A  is  the  maximum  dimension  over  all  its 
simplices,  or  is  infinite  if  the  maximum  does  not  exist.  The  vertex  set  V  of  A  is  the 
union  of  the  one-point  elements  of  A.  Fig.  4  shows  an  example  of  a  2-dimensional 


4  Notice  that  since  A  is  its  own  nonempty  subset,  the  simplex  A  is  also  a  face  of  A. 
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simplicial  complex.  A  subcollection  of  A  that  is  itself  a  simplicial  complex  is  called 
a  subcomplex  of  A.  A  subcomplex  of  A  is  the  p-skeleton  of  A,  denoted  by  Aip> ,  if 
it  is  the  collection  of  all  simplices  of  A  with  dimension  no  larger  than  p.  Thus,  the 
1-skeleton  is  the  underlying  graph  of  A. 

A  facet  of  a  simplicial  complex  A  is  a  maximal  face  of  A,  i.e.,  it  is  not  a  subset  of 
any  other  face.  A  simplicial  complex  is  connected  if  its  1-skeleton  (i.e.,  the  underlying 
graph)  is  connected  in  the  graph  sense. 

A  weighted  simplicial  complex  (WSC)  A  is  a  triple  ( V,S,w )5,  where  V  is  the 
set  of  vertices,  S  the  set  of  faces  of  A,  and  w  :  S  —>  {R+  U  {0}}  a  nonnegative 
weight  function  defined  for  each  face  in  S  with  w(v)  =  0  for  all  v  £  V.  We  define 
the  facet-only  weight  W f(A)  of  a  WSC  A  as 

Wf(A)  =  E  w&). 

RiE {facet  of  A} 


3  Minimum  Connected  Spanning  Subcomplex 

In  this  section,  we  show  that  the  MCSSub  problem  is  NP-complete,  and  we  pro¬ 
pose  two  approximation  algorithms  based  on  connected  set  cover  and  node-weighted 
Steiner  tree.  We  also  establish  the  approximation  ratios  of  the  two  algorithms  and 
analyze  their  time  complexity. 


3.1  NP-Completeness 

The  decision  version  (D-MCSSub)  of  the  MCSSub  problem  is  stated  as  follows:  let 
V(A)  denote  the  vertex  set  of  a  WSC  A  and  WF(A)  the  facet-only  weight  of  A. 
Given  a  WSC  A  =  (  V.  S.  w)  and  K  >  0,  is  there  a  connected  subcomplex  Asub 
of  A  such  that  V(Asub)  =  V  and  WF(Asub)  <  I\1  Then  we  have  the  following 
theorem. 

Theorem  1  The  D-MCSSub  problem  is  NP-complete. 

To  prove  the  NP-completeness,  we  reduce  a  classic  NP-complete  problem  -  the 
unweighted  set  cover  problem  to  the  MCSSub  problem. 

Proof  To  check  a  solution  to  the  D-MCSSub  problem,  we  only  need  to  verify  the 
following  points:  (i)  compute  the  facet-only  weight  of  the  solution  and  compare  the 
weight  with  K ;  (ii)  check  whether  the  underlying  graph  of  the  solution  is  connected  or 
not;  (iii)  check  whether  all  the  vertices  of  the  original  simplicial  complex  are  covered 
by  the  solution.  Since  all  these  can  be  done  within  polynomial  time,  the  D-MCSSub 
problem  is  NP. 

Given  is  an  unweighted  Set  Cover  instance  I,  i.e.,  a  universe  of  elements  U  and 
a  family  of  subsets  F  of  U.  For  each  element  u  G  U,  we  introduce  a  corresponding 
vertex  in  MCSSub  instance  I'.  We  introduce  one  additional  vertex  d.  For  each  set  /  G 


5  (<S,  w)  suffices  to  denote  the  WSC  since  V  C  <5.  but  we  use  the  redundant  ( V ,  S,  w)  for  convenience. 
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F,  we  introduce  a  corresponding  face  f  =  f  U  {d},  of  weight  1.  Being  a  simplcial 
complex,  all  subsets  of  f  are  also  introduced,  all  also  of  weight  1.  In  addition  to 
covering  all  vertices,  a  solution  to  /'  must  be  connected. 

Given  any  solution  SOL  to  I  of  cost  c,  we  can  construct  a  solution  to  I'  also  of 
cost  c.  For  every  set  /  in  SOL  not  including  d,  replace  /  with  /  U  {d}. 

On  the  other  hand,  given  any  solution  SOL'  to  I'  of  cost  c,  a  solution  to  I  of  the 
same  cost  can  be  constructed  as  follows:  for  any  face  /  £  SOL'  such  that  /  does  not 
appear  as  a  set  in  F,  replace  /  with  any  superset  of  /  —  {d}  appearing  in  F.  Note 
that  there  must  exist  at  least  one  such  superset. 

This  proof  is  for  general  D-MCSSub  problems.  It  can  be  shown  that  even  if  the 
weight  function  of  the  WSC  is  monotone  or  strictly  monotone6,  the  D-MCSSub  prob¬ 
lem  is  NP-hard.  But  it  is  still  possible  that  the  D-MCSSub  problem  under  some  spe¬ 
cial  structured  weight  function  is  P. 

In  the  following,  we  present  two  approximation  algorithms  for  the  MCSSub  prob¬ 
lem  both  with  performance  guarantee  0(ln  n),  where  n  is  the  number  of  vertices  in 
the  WSC.  Since  the  best  possible  approximation  ratio  for  the  set  cover  problem  is 
Inn  [5],  these  two  algorithms  are  order-optimal. 


3.2  Algorithm  Based  on  Connected  Set  Cover 

Let  A  be  a  set  with  finite  number  of  elements,  and  B  =  {Bi  C  A  :  i  =  1 , ....  n}  a 
collection  of  subsets  of  A  where  each  Bi  is  associated  with  a  weight  w(Bi)  >  0.  Let 
G  be  a  connected  graph  with  the  vertex  set  B.  A  connected  set  cover  (CSC)  Sc  with 
respect  to  (A,  B,  w ,  G)  is  a  set  cover  of  A  such  that  Sc  induces  a  connected  subgraph 
of  G.  The  minimum  connected  set  cover  (MCSC)  problem  is  to  find  the  CSC  with 
the  minimum  weight,  where  the  weight  of  a  CSC  Sc  is  defined  as 

w{Sc)  =  ^2  w{Bi). 

Bi^Sc 

From  a  WSC  A  =  (  V,  S.  w),  we  derive  an  auxiliary  undirected  graph  Ga  in  the 
following  way:  let  S  \  V  be  the  vertex  set  of  Ga,  and  connect  two  vertices  (non¬ 
vertex  faces  in  A)  Si  and  S2  if  and  only  if  Si  C  S2  f-  0  (i.e.,  Si  and  S2  have  at  least 
one  element  of  V  in  common).  Then  we  have  the  following  theorem  on  the  relation 
between  the  MCSSub  problem  and  the  MCSC  problem. 

Theorem  2  Let  A*  be  the  MCSSub  of  a  WSC  A  =  ( V.  S,  w)  and  Sf  the  MCSC  of 
(V,  S\V,w ,  Ga)-  Then  we  have 

wF(A*)  =  w(S£). 

Proof  The  proof  is  based  on  the  following  lemma. 

Lemma  1  Let  Sf  be  the  MCSC  of  ( V ,  S\V.  w.  Ga)-  For  any  face  S  £  Sf  with 
w(S)  >  0,  we  have  that  there  does  not  exist  a  face  S'  £  Sf,  such  that  S  C  S'. 


6  We  say  that  the  weight  function  satisfies  the  monotone  property  if  for  any  two  faces  S 1  C  S2- 
w(S  1)  <  w(S 2),  i  e„  the  weight  is  monotone  non-decreasing  with  respect  to  the  dimension  of  the  face. 
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Proof  (Proof  of  Lemma  1)  Suppose  that  for  some  face  S  £  Sf  with  w(S)  >  0, 
3  S'  £  < Sq  such  that  S  C  S'.  Let  S'c  =  Sf  \  s.  Obviously,  S'c  is  a  set  cover,  and 
w(S'c)  =  w(Sq)  —  w(S)  <  w(Sq).  On  the  other  hand,  since  S  O  S"  ^  0  implies 
S'  n  S'"  ^  0  for  any  face  S"  £  Sf,  it  follows  from  the  connection  rule  of  the 
auxiliary  graph  Ga  that  any  path  via  S  has  an  alternative  path  via  S'.  Thus,  S'c  is  a 
CSC,  leading  to  a  contradiction. 

Given  the  MCSC  Sf  of  ( V.  S  \  V.  w,  Ga),  we  can  obtain  a  connected  spanning 
subcomplex  A*c  by  mapping  each  element  of  Sf  to  a  face  in  A.  Since  the  facet-only 
weight  wf{A*c)  of  A*c  only  counts  facets  in  A*c,  it  follows  that  wf{Aq)  <  w(Sf). 
Based  on  Lemma  1,  we  have  that  every  element  of  Sf  with  positive  weight  is  a  facet 
in  A*c,  and  thus 

wf(A*)  <  wf(A*c)  =  w(Sf), 

where  A*  is  the  MCSSub  of  A. 

On  the  other  hand,  the  facets  of  A*  leads  to  an  CSC  Sf,  and  w(Sf)  =  wf(A*). 
It  implies  that 

w(Sf)<w(S*A)  =  wF(A*). 

Thus,  wf(A*)  =  w(Sf). 

3.2.1  Algorithm 

Based  on  Theorem  2,  we  can  reduce  the  MCSSub  problem  of  a  WSC  A  =  (V,  S,  w ) 
to  the  MCSC  problem  (V,  S  \V,w,  Ga)-  We  obtain  the  following  Set  Cover  based 
Algorithm  (SCA)  for  the  MCSSub  problem. 


Algorithm  1  SCA  for  MCSSub: 

INPUT:  A  WSC  A  =  ( V. ,  S,  w). 

OUTPUT:  An  approximate  MCSSub  Ac  of  A. 

1.  Derive  the  auxiliary  graph  Ga- 

2.  Find  an  approximate  MCSC  Sc  of  (V,  S  \V,w,  Ga)  by  using  the  greedy  algo¬ 
rithm  for  MCSC  (Algorithm  2). 

3.  Transform  Sc  to  a  connected  spanning  subcomplex  Ac  by  mapping  each  element 
of  Sc  to  a  face  in  A. 


Zhang  et  al.  propose  a  greedy  approximation  algorithm  for  the  unweighted  MCSC 
problem  [23],  i.e.,  w{Bf)  =  1  for  all  i.  The  original  algorithm  in  [23]  has  a  flaw  and 
the  established  approximation  ratio  is  incorrect.  In  [20],  the  flaw  is  corrected  and 
a  stronger  result  on  the  approximation  ratio  is  shown.  By  generalizing  their  greedy 
approach,  we  develop  a  greedy  algorithm  for  the  weighted  MCSC  problem. 

Before  stating  the  algorithm,  we  introduce  the  following  notations  and  definitions. 
For  two  sets  Si,  62  £  S,  let  dist(;(,5'i ,  S2)  be  the  length  of  the  shortest  path  between 
Si  and  S2  in  an  auxiliary  graph  G,  where  the  length  of  a  path  is  given  by  the  number 
of  edges;  Si  and  S2  are  said  to  be  graph-adjacent  if  they  are  connected  via  an  edge 
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in  G  (i.e.,  distG(Si,  S2)  =  1),  and  they  are  said  to  be  cover-adjacent  if  Si  C\S2  ^  0. 
Notice  that  in  a  general  MCSC  problem,  there  is  no  connection  between  these  two 
types  of  adjacency.  The  cover-diameter  Dq(G)  is  defined  as  the  maximum  distance 
between  any  two  cover-adjacent  sets,  i.e., 

DC(G)  =  max{distG(S'i,  S2)  \  SltS2  £  S  and  Si  n  S2  ±  0}. 

For  the  MCSC  problem  derived  from  the  MCSSub  problem  of  a  WSC  A,  we  have 
that  Dc(Ga)  =  1- 

At  each  step  of  the  algorithm,  let  TZ  denote  the  collection  of  the  subsets  (faces 
of  A)  that  have  been  selected,  and  U  the  vertex  subset  of  A  that  has  been  covered. 
Given  TZ  ^  0  and  a  set  S  £  S  \  TZ,  an  1Z  — >  S  path  is  a  path  {So,  Si, ...,  Sk}  in  G 
such  that  (i)  So  £  TZ',  (ii)  Sk  =  S;  (iii)  Si,...,  Sk  £  S\TZ.  We  define  the  weight  ratio 
r(Ps)  of  Ps  as 

tp  \  __  w(S(Ps)\TZ)  _  YlseS(ps)\nw(s ) 

1  \VN(Ps)\  IMPs)|  ’ 

where  S(Ps)  \  TZis  the  subsets  (faces  in  S)  of  Ps  that  are  not  in  7 Z,  and  |Vjv(-Pg)|  is 
the  number  of  vertices  of  A  that  are  covered  by  Ps  but  not  covered  by  7 Z. 


Algorithm  2  A  Greedy  Algorithm  for  MCSC. 

INPUT:  ( V,S\V,w,Ga ) 

OUTPUT:  A  CSC  K. 

1.  Choose  So  £  S\V  such  that  the  weight  ratio  r  (So)  defined  in  (1)  is  the  minimum, 
and  let  TZ  =  {So}  and  U  =  So- 

2.  WHILE  V  \U  ^0  DO 

2.1.  For  each  S  £  S  \  (V  U  TZ)  which  is  cover-adjacent  or  graph-adjacent  with  a 
set  in  TZ,  find  a  shortest 1  TZ  — >  S  path  Ps- 

2.2.  Select  Ps  with  the  minimum  weight  ratio  r(Ps)  defined  in  (1),  and  let  TZ  = 
TZ  U  Ps  (add  all  the  subsets  of  Ps  to  TZ)  and  U  =  U  U  Vn(Ps)- 

END  WHILE 

3.  RETURN  TZ. 


3.2.2  Approximation  Ratio 

The  approximation  ratio  of  SCA  is  determined  by  Step  2,  i.e.,  the  approximation  ratio 
of  the  greedy  algorithm  for  the  MCSC  problem.  First,  we  establish  the  following 
lemma. 

7  Notice  that  the  shortest  path  is  defined  in  terms  of  the  number  of  edges,  not  the  total  weight  of  all 
vertices  along  the  path. 
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Lemma  2  Given  a  weighted  MCSC  problem  (V.  S  \V,w,  G)  with  Dc(G)  =  1,  let 

R w  =  s*s,  ,qn.  (2) 
min  w(b )  1 
SgS 

Then  the  approximation  ratio  of  the  greedy  algorithm  for  MCSC  is  at  most  Rw  + 
H(  7  —  1),  where  7  =  max{|S'|  |  S  £  S  \  V}  is  the  maximum  size  of  the  subsets  in  S 
and  Hf)  is  the  harmonic  function. 

Proof  (Proof  of  Lemma  2)  The  proof  is  based  on  the  classic  charge  argument.  Let  S* 
be  an  optimal  solution  to  the  weighted  set  cover  problem  [V.  S  \  V.  w),  and  TZ  the 
solution  returned  by  the  greedy  algorithm  for  the  weighted  MCSC  problem  (V,  S  \ 
V,w,G)  with  Dc(G)  =  1.  Let  w(S*)  and  w(JZ)  denote  the  total  weight  of  the 
subsets  included  in  S*  and  TZ,  respectively.  In  the  following,  we  will  show  that 

+  ,3) 


Let  Sf  be  an  optimal  solution  to  the  weighted  MCSC  problem  ( V.  S  \V,w,  G ) 
with  Dc(G)  =  1.  Since  w(S*)  <  w(Sf),  Lemma  2  follows  immediately  from  (3). 
To  prove  (3),  we  apply  the  classic  charge  argument.  Each  time  a  subset  So  (at  step 

1)  or  a  shortest  TZ  — >  S  path  Pg  (at  step  2)  is  selected  to  be  added  to  7 Z,  we  charge 

each  of  the  newly  covered  elements  (at  step  1)  or  r(P|)  defined  in  (1)  (at  step 

2) .  Notice  that  when  Dc{G)  =  1,  the  shortest  TZ  S  path  Pg  is  only  a  single  edge 
connecting  some  subset  in  TZ  and  S,  and 

rp*3  =  ™{S{Ps  \  n))  =  w(s) 

1  s>  |Viv(P|)|  \s\uy 

During  the  entire  procedure,  each  element  of  V  is  charged  exactly  once.  Assume 
that  step  2  is  completed  in  K  —  1  iterations.  Let  Pf  be  the  shortest  TZ  — >•  S  path 
selected  by  the  algorithm  at  iteration  i.  Let  C(v)  denote  the  charge  of  an  element  v 
in  V.  Then  we  have  that 


K- 1 

=  E 

E 

C(v) 

vGV 

2=0  ^ 

»e  vn(p 

K- 1 

=  E 

2=0  1 

E 

veSi\u 

w(Si) 
|5i  \  U | 

K- 1 

=  E 

2=0 

w(Si)  -- 

=  w(JZ), 

where  P*Q  =  {SQ}. 
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Suppose  that  S*  =  {S^, is  a  minimum  weighted  set  cover  for  [V.  S  \ 
V,w}.  Since  an  element  of  V  may  be  contained  in  more  than  one  subset  of  S*,  it 
follows  that 

N 

E  c («)  <  E  E  c («)•  (5) 

v(zV  i— 1  v£S* 

Next  we  will  show  an  inequality  which  bounds  from  above  the  total  charge  of  a 
subset  in  S*,  i.e.,  for  any  B*  e  S*, 

E  °(v)  <  lR™  +  H(\s*\  -  VMS*).  (6) 

v£B* 

Let  m  (i  =  0, 1, K)  be  the  number  of  elements  of  S*  that  have  not  been 
covered  by  S  after  iteration  i  —  1,  where  step  1  is  considered  as  iteration  0.  Let 
{ii, ik}  denote  the  subsequence  of  {i  =  0, 1, K  —  1}  such  that  rii  —  rii+i  >  0. 
For  each  element  a  covered  at  iteration  *i,  if  =0,  based  on  the  greedy  rule  at  step 
1,  we  have  that 

C\v)  =  r{P*So)  <  ^ P (7) 

i 

Otherwise, 


C(v)  =  r(P*Sii) 


wjS^)  <  w(S*)Rw 

|5ii  \  U | 


(8) 


The  inequality  in  (8)  is  due  to  the  fact  that  S,n  covers  at  least  —  7i(il+1)  elements 
of  V,  i.e.,  |5ij  \  U\  >  riij  —  Summing  up  (7)  and  (8), 


C(v)< 


w(S*)Rw 

ni± 


(9) 


Consider  two  cases: 

(i)  If  all  the  elements  of  S*  have  been  covered  after  iteration  i\,  i.e.,  n^1+ 1)  =  0, 
then 

E  C(v)  ^  E  W{S*)Rw  =  w(S*)Rw.  (10) 

4^  4£  no 

(ii)  If  not  all  the  elements  of  S*  have  been  covered  by  TZ  after  iteration  i i,  S*  becomes 
cover-adjacent  with  1Z  and  thus  a  candidate  for  being  selected  at  the  following 
iterations.  At  each  iteration,  for  each  element  v  €  S*  covered  at  iteration  ij 
(j  =  2, k),  the  greedy  rule  at  step  2  still  yields 

C(v)=r(P*Si_)<r(P*s,) 

3 

w(S*)  w{S*) 


(ii) 
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It  follows  from  (9,1 1)  that 


E  C(v'>  -  w(s*)(nh 

ves * 


Tl(*i+1)) 

IL' 


1 


k  1 

+w(S*)Y^ K  -«(<j+i))  — 


3= 2 


=  w(S*) 


(12) 


Here  we  have  used  the  fact  that  =  nj.  .+1)  •  It  is  because  between  iteration 

ij  and  iteration  no  elements  of  S*  are  covered. 

For  the  summation  term  in  (12),  we  have  the  following  inequality: 


£: 

3= 2 


^E 


i 


rij . 

i=2 


nHi+i)  +  1 


=  H(nl2)<H(\S*\-l). 

The  last  inequality  is  due  to  the  fact  that  rii2  <  —  1  =  |5*|  —  1. 

Eqn.  (6)  is  a  direct  consequence  of  (10),  (12),  and  (13).  Thus,  using  (4-6), 


N 


w{k)  =  E  -EE  c(v) 


vev 

N 


i=  1  v£S? 


2=1 


<  [Rw  +  H(p/  —  l)]tu(iS*). 


03) 


Then,  as  a  direct  consequence  of  Lemma  2,  we  have  the  following  theorem  on  the 
approximation  ratio8  of  the  greedy  algorithm  for  the  MCSC  problem  with  Dc{G)  = 

1. 


Theorem  3  Let  A*  be  the  MCSSub  of  a  WSC  A  =  ( V.  S.  w)  and  Ac  be  the  solution 
returned  by  Algorithm  1.  Let  Rw  be  defined  as  in  (2).  Then  we  have 

<  Rw  +  H(dimA), 

where  dim  A  is  the  dimension  of  A  and  Hf)  is  the  harmonic  function. 

From  Theorem  3,  we  see  that  the  approximation  ratio  depends  on  the  ratio  Rw  of 
the  maximum  weight  to  the  minimum  weight.  It  is  shown  in  the  following  theorem 
that  if  Rw  is  unbounded,  then  the  scaling  order  of  the  approximation  ratio  can  be  as 
bad  as  linear  with  respect  to  the  number  of  vertices  in  the  simplicial  complex. 

8  The  approximation  ratio  of  the  greedy  algorithm  for  general  weighted  MCSC  problem  is  still  an  open 
problem. 
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Theorem  4  Let  n  be  the  number  of  the  vertices  in  a  WSC  A  =  (V,  S,  w),  and  Rw 
defined  as  in  (2).  If  Rw  is  unbounded ,  then  the  approximation  ratio  of  Algorithm  1 
for  the  MCSSub  problem  of  A  is  f2(n). 

Proof  Consider  a  specific  example:  A  is  a  (n  —  l)-dimensional  simplex  with  the 
vertex  set  V  =  {vi . vn},  and  all  the  weights  of  the  faces  are  infinite  except  for  the 
following  five  faces: 


w(Si)  =  W  ({til,..., un  })  = 

w(S2)  =U>({ui,».,Uf,U(»+2)})  =  2  ’ 
w(S3)=w({v{i+l),...,v{9+l)})=\, 
w(S4)  =  W  ({un,...,u„})  = 
w(S5)  =  W  ({u(„+1),...,Un})  =  1- 

For  ease  of  presentation,  we  have  assumed  that  n  is  a  multiple  of  4.  By  applying 
Algorithm  1,  we  reduce  the  MCSSub  problem  for  A  to  the  MCSC  problem  (  V,  S  \ 
V,  w,  Ga )■  Due  to  the  weight  assignment,  it  suffices  to  only  consider  the  subgraph  of 
G a  induced  by  the  above  five  faces,  as  shown  in  Fig  5. 

The  optimal  solution  A*  to  the  MCSSub  problem  is  given  by 

A*  =  {S  g  5  |  S  C  ^  or  S3  or  S5}, 


and 

wf(A*)  =  w(S'2)  +  w(S3)  +  w(S5)  =  2. 

On  the  other  hand,  the  solution  Ac  returned  by  Algorithm  1  is  given  by 


Ac  =  {SeS\SCS1  or  S4}, 


and 

77 

wf{Ac)  =  w(Si)  +  w(S4)  =  -  +  -■ 

Z  o 

Specifically,  Si  is  firstly  selected,  and  then  ,S'i .  Thus, 

wf{Ac)  n  1 

wf(A*)  16  4  y  J 

It  follows  that  the  approximation  ratio  of  Algorithm  1  is  Q{n). 


From  Theorem  4,  we  see  that  Algorithm  1  is  not  suitable  for  the  MCSSub  problem 
of  a  WSC  A  if  its  weight  function  has  a  relatively  wide  range.  As  shown  next  in 
Sec.  3.3,  the  other  approximation  algorithm  based  on  the  Steiner  tree  does  not  have 
this  issue:  its  approximation  ratio  does  not  depend  on  the  range  of  the  weight  function. 
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3.3  Algorithm  Based  on  Steiner  Tree 

From  a  WSC  A  =  (V.  S.  w ),  we  derive  an  undirected  graph  H a  with  the  vertex  set 
S:  for  each  face  S  €  S\V  (i.e.,  the  faces  that  are  not  the  vertices  of  A),  we  replace 
it  by  a  vertex  v s  in  Ha  and  connect  vs  to  all  the  vertices  of  S.  The  weight  w(vs) 
assigned  to  the  vertex  vs  is  the  weight  w(S)  of  the  face  S.  Notice  that  the  weight 
of  vertices  in  Ha  corresponding  to  the  vertices  in  A  {i.e.,  V)  is  zero.  Fig.  6  shows 
an  example  of  the  derivation  of  the  graph  from  a  2-simplex.  We  have  the  following 
theorem  on  the  relation  between  the  MCSSub  of  A,  the  Steiner  tree  of  II a  that  spans 
the  vertex  set  V  of  A  and  the  minimum  connected  dominating  set9  of  II a- 

Theorem  5  Let  A*  denote  the  MCSSub  of  a  WSC  A  =  ( V. ,  S,  w),  T*  the  Steiner  tree 
of  Ha  that  spans  the  vertex  set  V  of  A,  and  Df,  the  minimum  connected  dominating 
set  of  Ha-  Then  we  have  that 

wF{A*)=w{T*)=w{D*c). 

Proof  First  we  show  that  wF{A*)  =  w(T*).  Since  every  connected  spanning  sub¬ 
complex  A'  of  A  corresponds  to  a  connected  subgraph  of  II a  which  only  contain- 
s  the  vertices  of  A  and  the  vertices  representing  the  facets  of  A',  it  follows  that 
w(T*)  <  wf(A*).  On  the  other  hand,  since  by  contradiction,  there  is  a  one-to-one 
mapping  between  the  vertices  of  the  Steiner  tree  of  Ha  and  the  vertices  plus  the 
facets  of  a  connected  spanning  subcomplex  of  A,  it  follows  that  wF{A *)  <  w{T*). 

Next  we  show  that  u>(T*)  =  w{Df).  Notice  that  the  vertex  set  V  of  A  is  a 
dominating  set  of  Ha-  Since  the  Steiner  tree  T*  of  Ha  spans  the  vertex  set  V,  T*  is 
a  CDS  of  Ha-  Thus,  w{Dq)  <  w(T*).  On  the  other  hand,  given  the  minimum  CDS 
I)*c  of  Ha,  since  each  vertex  v  in  the  vertex  set  V  is  either  in  !)},  or  a  neighbor  of 
some  face  in  I)  f  and  the  weights  of  the  vertices  in  V  are  all  zero,  the  combination  of 
V  and  Dq  yields  a  connected  subgraph  of  Ha  that  spans  V  with  the  same  weight  as 
D*c.  Thus,  w (T*)  <  w{D*c). 

Based  on  Theorem  5,  we  propose  the  following  Steiner  Tree  based  Algorithm 
(STA)  for  the  MCSSub  problem. 


Algorithm  3  STA  for  MCSSub: 

INPUT:  A  WSC  A  =  {V,  S,  w). 

OUTPUT:  An  approximate  MCSSub  Ac  of  A. 

1.  Derive  the  graph  Ha  from  A. 

2.  Obtain  an  approximate  Steiner  tree  T  of  Ha  by  using  the  algorithms  given  in  [7, 

8]. 

3.  Transform  T  to  a  connected  spanning  subcomplex  Aq  of  A  by  mapping  each 
vertex  ofT  to  a  face  of  A. 

9  A  dominating  set  of  a  graph  is  a  subset  of  vertices  such  that  every  vertex  of  the  graph  is  either  in  the 
subset  or  a  neighbor  of  some  vertex  in  the  subset,  and  a  connected  dominating  set  (CDS)  is  a  dominating 
set  where  the  subgraph  induced  by  the  vertices  in  the  dominating  set  is  connected.  The  CDS  problem  asks 
for  a  CDS  with  the  minimum  total  weight,  and  it  is  shown  to  be  a  special  case  of  the  MCSC  problem  [20], 
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Since  approximation  only  occurs  in  Step  2,  the  approximation  ratio  of  STA  is 
equal  to  that  of  the  algorithm  for  the  node-weighted  Steiner  tree  problem.  The  best 
approximation  ratio  is  known  to  be  (1.35  +  e)  Inn  for  any  constant  e  >  0,  where  n 
is  the  number  of  vertices  of  A  and  is  also  the  number  of  terminals  in  the  Steiner  tree 
of  Ha  [7].  Here  we  do  not  try  to  find  the  CDS  I)*c  of  Ha  at  step  2,  because  the  best 
known  approximation  ratio  for  the  CDS  problem  is  (1.35  +  e)  In  n( /lAj  [6,  7].  Since 
rii  a A  j  3>  n,  the  latter  approximation  ratio  is  much  worse  than  the  former  one. 


3.4  Time  Complexity  Analysis 

Here  we  analyze  the  time  complexity  of  SCA  and  STA  for  the  MCSSub  problem. 
Given  a  WSC  A  =  (V,  S,  w ),  let  n  =  \V\  denote  the  number  of  vertices  in  A,  m  = 
|<S  \  V |  the  number  of  non-vertex  faces  in  A,  and  d  the  dimension  of  A.  Recall  that 
the  existence  of  edges  in  the  auxiliary  graph  Ga  for  SCA  and  the  derived  graph  II a 
for  STA  depends  entirely  on  whether  the  two  non- vertex  faces  overlap  and  whether 
the  vertex  is  contained  in  the  non-vertex  face,  respectively.  It  implies  that  all  the 
information  of  these  two  graphs  can  be  easily  retrieved  from  the  WSC  A.  Thus,  Step 
1  in  both  algorithms  can  be  skipped  in  the  implementation,  and  the  time  complexity 
of  both  algorithms  is  determined  by  their  Step  2. 

Step  2  of  SCA  is  to  apply  the  greedy  algorithm  to  the  MCSC  problem  (V,  S  \ 
V ,  w,  Ga)-  It  takes  0(m)  time  to  complete  Step  1  of  the  greedy  algorithm.  Since  at 
least  one  vertex  becomes  covered  at  each  iteration  of  Step  2  of  the  greedy  algorithm, 
there  are  at  most  n  —  1  iterations.  At  each  iteration,  the  weight  ratios  of  at  most  to 
faces  are  computed,  and  due  to  the  fact  that  cover-adjacent  faces  are  graph-adjacent, 
the  weight  ratio  of  each  face  is  done  in  constant  time.  Thus,  the  running  time  of  SCA 
is  0(m  +  nm )  =  0(nm). 

Since  the  derived  graph  Ha  has  n  +  m  vertices  and  0(dm )  edges  and  the  Steiner 
tree  has  n  terminals  to  cover,  it  follows  from  [13]  that  the  running  time  of  Step  2  of 
STA  is  0(dnm2  +  nm2  log  to).  From  the  above,  we  see  that  the  time  complexity 
of  STA  is  significantly  higher  than  that  of  SCA.  This  is  mostly  because  the  approx¬ 
imation  algorithm  for  the  Steiner  tree  requires  the  computation  of  the  shortest  paths 
between  all  vertex  pairs. 

We  point  out  that  while  the  Steiner  tree  based  algorithm  has  a  higher  complexity, 
it  can  offer  better  performance  in  a  WSC  with  a  large  weight  range.  In  a  simulation 
example  of  random  simple  complexes,  we  consider  a  case  where  each  face  weight 
takes  only  two  values  wmin  and  wmax  with  equal  probability.  With  wmin  =  1, 
w-max  =  10000,  and  1000  Monte  Carlo  runs  for  a  200- vertex  random  simplicial 
complex10  [14],  we  find  that  the  total  weight  of  the  solution  returned  by  the  set  cov- 

10  A  random  simplicial  complex  A(n,D,  p)  with  n  vertices,  dimension  at  most  V).  and  a  D- 
dimensional  probability  vector  p  =  {pi , P2 ,  ■  ■  ■ ,  PD  }  is  constructed  in  a  bottom-up  manner:  first  n  vertices 
are  fixed,  which  are  the  0-simplices  of  A ,  and  then  higher-dimensional  simplices  are  generated  inductively. 
Specifically,  for  each  1  <  i  <  D,  after  all  the  simplices  with  dimension  lower  than  i  have  been  gener¬ 
ated,  consider  every  i-tuple  of  vertices:  if  they  have  formed  all  the  lower  dimensional  simplices,  then  an 
/-simplex  consisting  of  them  is  generated  with  probability  pi.  Notice  that  a  random  simplicial  complex 
A  (n,  1,  p)  is  the  random  graph  introduced  by  Erdos  and  Renyi  [4], 
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er  based  algorithm  can  be  1.7  times  that  of  the  solution  returned  by  the  Steiner  tree 
based  algorithm.  These  two  algorithms  thus  offer  a  tradeoff  between  performance  and 
complexity. 


4  Simulation  Results 


In  this  section,  we  present  simulation  results  on  the  performance  of  the  two  approxi¬ 
mation  algorithms  (SCA  and  STA)  for  the  broadcast  problem  in  an  MR-MC  network. 
We  consider  a  dense  MR-MC  network,  where  all  the  nodes  are  within  each  other’s 
transmission  range,  and  we  aim  to  minimize  the  total  energy  consumption  of  the 
broadcast. 

There  are  12  non-overlapping  channels  /,;  (1  <  i  <  12),  possibly  with  differ¬ 
ent  communication  rates  r,,  available  for  the  MR-MC  network,  and  each  node  is 
equipped  with  4  radios.  At  the  beginning  of  the  broadcast,  each  node  randomly  se¬ 
lects  4  of  the  12  channels  for  its  4  radios.  As  discussed  in  Sec.  1.3,  the  nodes  which 
share  at  least  one  common  channel  form  a  clique,  and  there  is  a  one-to-one  corre¬ 
spondence  between  the  cliques  and  the  faces  of  the  derived  WSC.  The  weight  of  the 
face  is  defined  as  the  energy  consumption  of  the  broadcast  within  the  corresponding 
clique,  i.e.,  the  sum  of  the  transmission  energy  and  the  reception  energy.  Let  S  be  a 
face  containing  k  +  1  nodes  and  {fsj  :  j  =  1,  2, ...,  g}  the  q  (1  <  q  <  12)  common 
channels  shared  by  the  k  +  1  nodes.  Assume  that  if  a  node  in  the  clique  is  selected  as 
relay,  it  will  choose  the  common  channel  with  the  maximum  communication  rate  to 
transmit.  Then  the  weight  w(S )  of  the  face  S  is  given  by 


w(S) 


( Ptx  +  kPrx) 


L 

max  { rSj }’ 
.7=1, •••>9 


where  Ptx  and  Prx  are  the  transmission  power  and  the  reception  power,  respectively, 
and  L  is  a  constant. 

In  Fig.  7,  the  average  total  energy  of  the  solutions  returned  by  SCA  and  STA  is 
compared  with  that  of  the  MST  with  respect  to  the  underlying  graph  of  the  WSC. 
The  average  is  taken  over  10  random  channel  assignments.  Notice  that  although  two 
different  links  on  the  same  channel  are  treated  separately  when  the  MST  is  derived, 
the  transmission  energy  corresponding  to  them  is  counted  only  once  to  exploit  the 
wireless  broadcast  advantage  when  the  total  energy  of  the  MST  is  computed.  We  see 
that  the  performances  of  SCA  and  STA  are  extremely  close,  and  their  performances 
are  significantly  better  than  that  of  MST. 


5  Conclusion  and  Future  Work 

In  this  paper,  we  study  the  minimum  cost  broadcast  problem  in  multi-radio  multi¬ 
channel  ad  hoc  networks,  where  the  total  cost  is  the  sum  of  the  costs  associated  with 
the  transmissions  during  the  broadcast.  We  formulate  it  as  the  minimum  spanning 
problem  in  simplicial  complexes.  We  show  that  it  is  NP-complete.  Hence  we  pro¬ 
pose  two  approximation  algorithms  for  this  minimum  spanning  problem:  one  is  to 
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transform  it  into  the  connected  set  cover  problem;  the  other  is  to  transform  it  into 
the  node-weighted  Steiner  tree  problem  and  then  apply  the  corresponding  algorith- 
m.  Despite  their  distinct  approaches,  both  approximation  algorithms  are  shown  to  be 
order-optimal  and  offer  a  tradeoff  in  terms  of  performance  vs.  complexity. 

As  a  starting  point,  we  have  assumed  that  the  channel  assignment  scheme  is 
designed  independent  of  the  broadcast  scheme.  The  joint  optimization  of  the  two 
schemes  will  further  reduce  the  broadcast  cost.  Another  future  direction  is  to  develop 
distributed  versions  of  the  approximation  algorithms  for  the  minimum  cost  broadcast 
problem. 
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V\  V2  VI  V2 

(a)  Graph  (b)  Simplicial  Complex 

Fig.  1  Graph  and  simplicial  complex:  V  =  {uo,  vu  ^2}*  £(a)  —  {(^0 ,  vi),  (vo,  V2),  (v±  ,  ^2)}*  S(b)  — 
{(no,  vi,  V2),  (no,  ni),  (no,n2),  (ni,n2),  (no),  (ni),  (n2)}- 


Fig.  2  An  illustration  of  an  MR-MC  network  and  the  constructed  simplicial  complex.  The  parameters 
within  the  braces  are  the  channels  which  each  node  can  access.  In  the  communication  graph  derived  from 
the  network,  a  link  exists  between  two  nodes  if  and  only  if  two  nodes  are  within  each  other’s  transmission 
range  and  they  share  at  least  one  common  channel.  Notice  that  a  clique  in  the  communication  graph  may 
not  be  a  clique  in  the  MR-MC  network  (correspondingly,  a  simplex  in  the  simplicial  complex),  e.g.,  the 
three  nodes  of  the  right  empty  triangle  (they  do  not  share  a  common  channel). 


20 


W.  Ren,  Q.  Zhao,  R.  Ramanathan,  J.  Gao,  A.  Swami,  A.  Bar-Noy,  M.  Johnson,  R  Basu 


vo 


Fig.  3  A  simplicial  complex  whose  MCSSub  is  the  three  filled  triangles,  and  is  not  a  “tree”  (the  integers 
are  the  weights  of  the  simplices). 


V5 


Fig.  4  A  simplicial  complex  A 
plices: 

{«1 ,  V2}.{V2,V3},{V3,V4,},{V3,V5},{V4,V5}), 


with  6  vertices  (0-dimensional  sim- 
5  edges  (1 -dimensional  simplices: 

and  1  filled  triangle  (2-dimensional  simplex: 
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Si 


Fig.  5  The  subgraph  of  G  a  induced  by  the  five  faces  S\,  S2,  S3,  S4,  and  S3  with  finite  weights. 


zA  Ha 

Fig.  6  The  derived  graph  of  a  2-simplex  (squares  in  Ha  represent  the  faces  that  are  not  vertices  of  zA). 


Fig.  7  Average  total  energy  vs.  number  of  nodes.  Parameters:  Ptx=l,  Prx  =  0.01,  L  =  100,  r{  =  i  for 
1  <  i  <  12. 


